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Abstract 

The sensitivity of a search for sources of TeV neutrinos can be improved by grouping 
potential sources together into generic classes in a procedure that is known as source 
stacking. In this paper, we define catalogs of Active Galactic Nuclei (AGN) and use 
them to perform a source stacking analysis. The grouping of AGN into classes is 
done in two steps: first, AGN classes are defined, then, sources to be stacked are se- 
lected assuming that a potential neutrino flux is linearly correlated with the photon 
luminosity in a certain energy band (radio, IR, optical, keV, GeV, TeV). Lacking 
any secure detailed knowledge on neutrino production in AGN, this correlation is 
motivated by hadronic AGN models, as briefly reviewed in this paper. 

The source stacking search for neutrinos from generic AGN classes is illustrated 
using the data collected by the AMANDA-II high energy neutrino detector during 
the year 2000. No significant excess for any of the suggested groups was found. 
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1 Introduction 

With AMANDA-II, different searches for high energy neutrino point sources 
have been performed [1,2,3], showing that neither the all sky search for hot 
spots in the neutrino sky nor the search at predefined and astrophysically 
motivated source positions has led to a discovery of a statistically significant 
neutrino signal. This however docs not exclude the possibility that the sig- 
nal of a superposition of generically equal sources, each contributing below 
the individual significance threshold, sums up to a significant signal for that 
specific source type. In case of no detection, the stacked signal may be used 
to obtain a limit on neutrino emission from that generic source type. Source 
stacking methods have been applied in gamma-ray astronomy, e.g., in searches 
for non-blazar AGN classes at GeV photon energies [4,5]. 

AGN belong to the most promising potential neutrino sources. Since the de- 
tailed mechanism of neutrino production in AGN is still unknown, in this paper 
a neutrino-production-model free attempt to classify the sources according to 
their geometrical properties and electromagnetic emission is introduced. Based 
on this classification, we develop an AGN stacking analysis to search for TeV 
neutrinos. 

In Sec. 2, we describe in detail the adopted systematic classification of AGN 
based on an axisymmetric model [6]. Possible neutrino production scenarios 
in AGN arc reviewed in Sec. 3. Pion production and subsequent decay results 
in correlated neutrino and photon production. Neutrino fiuxcs from these re- 
actions are expected to be of the same order of magnitude as the gamma-ray 
fluxes. However, photons interact in the source environment possibly leading 
to a different photon spectrum than the neutrino spectrum. In Sec. 4, we de- 
fine catalogs of AGN classes based on our classification, and select from them 
interesting neutrino candidate sources. Hence, we order sources according to 
their photon fiux at different energies. In Sec. 5, a procedure is presented to 
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determine the optimum number of sources from each catalog to be included in 
the stacking analysis. This procedure is applied to determine source samples 
to be analyzed with the AMANDA neutrino telescope, located at the South 
Pole [7] . Thus selected sources in AMANDA's field of view are in the Northern 
sky. Assuming hypothetical distributions of the source strength as a function 
of the ordered source number we determine under which conditions a gain in 
sensitivity is obtained with the source stacking method with respect to other 
methods [1,8]. In Sec. 6, we illustrate our method by evaluating the signal 
from the resulting source samples using a data set collected by AMANDA in 
the year 2000 [1]. 



2 AGN classification 

In this section we describe a possible classification of AGN useful to define 
catalogs of sources to be stacked. Historically, a large variety of AGN (like 
Seyfert galaxies, radio galaxies or quasars) have been named due to their ap- 
pearance from Earth-based telescopes. The observational differences among 
the various AGN types can be partially explained with a geometrically ax- 
isymmetric model as the result of different inclination angles, defined as the 
angle between the line of sight and the AGN axis [6]. A pictorial scheme of 
an AGN showing the basic ingredients of the axisymmetric model is shown in 
Fig. 1. It consists of a rotating supermassive black hole, two jets with matter 
outflowing along the rotation axis and an accretion disk of matter perpendic- 
ular to the rotation axis. Aside from these geometrical differences, others will 
be outlined below. 

We describe the AGN classiflcation illustrated in Fig. 2 based on the mor- 
phology of the host galaxy, the luminosity and the inclination angle. AGN are 
generally divided into radio-loud and radio-weak sources as indicated in the 
first branching of the scheme. Radio-loud AGN (left branch of Fig. 2) can be 
classified according to radio jet lengths: compact objects, where the jets get 
stuck in dense matter [9], are discussed with more detail in Sec. 2.1 and AGN 
with fully evolved jets of a length of 100 kpc up to several Mpc in Sec. 2.2. 

Typically, AGN spectra are composed of a thermal part, the so-called blue 
bump with the maximum at optical-UV frequencies, and a non-thermal part 
extending over up to 20 orders of magnitude in frequency. The blue bump is 
interpreted as thermal radiation of the warm inner accretion disk. The low 
energy component of the non-thermal spectrum in the radio to soft X-ray 
regime is assumed to be due to synchrotron radiation of electrons gyrating in 
a magnetic field. The origin of the high-energy component, if present, can be 
explained by hadronic or leptonic models (see Sec. 3). The non-thermal high 
energy photon emission is known to be highly variable on a wide range of time 
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scales, from less than one hour to months [10] . Within the thermal spectrum, 
some emission lines can be found. 

Emission lines are classified as broad lines or as narrow lines. The axisymmetric 
model assumes that broad lines are due to fast moving dust clouds {v ~ 1000 
km/s— 5000 km/s) near the accretion disk. On the other hand, narrow lines 
originate in slowly moving dust clouds (y ~ 500 km/s) outside the torus. 
Depending on the angle of observation, the broad line emission from the clouds 
at the center may or may not be hidden by the torus. This model is confirmed 
by broad line observations in scattered light from some AGN, which show in 
direct light only narrow lines [11]. 

The radio emission of AGN is assumed to originate mostly in relativistic jets 
where it is caused by synchrotron radiation of electrons moving along the 
jet. AGN are called radio-loud if the ratio of the radio flux at 5 GHz to the 
optical flux is larger than 10 [6]. Radio-loud AGN are located in elliptical 
galaxies, while radio- weak AGN are located in spiral galaxies [12], and rarely 
in ellipticals. 

Observations and jet models [13] show a flat radio spectrum for the flux density 
F {F (X u", a > —0.5) ^ for the radio core, i.e. for the inner part of the jet. In 
contrast, radio lobes and hot spots located at the outer end of the jet typically 
show a steeper spectrum with spectral indices from —0.5 to —0.6 in hot spots 
and from —0.8 to —1.0 in lobes. 

For the radio- weak AGN (right branch of Fig. 2), a luminosity dependent 
classification in optical wavelengths yields the division into quasars and Seyfert 
galaxies. The intrinsically stronger objects are the radio-weak quasars and 
the Radio Intermediate Quasars (RIQ) [14]. The radio- weak quasars are seen 
from moderate inchnation angles (20° — 60°), the RIQ are interpreted as the 
same objects, but seen from smaller angles. The radio emission of RIQ is 
relativistically beamed, similarly to blazars (see Sec 2.2) [14]. The weaker 
objects are the Seyfert galaxies, classified as Seyfert 1 galaxies, if the core and 
the broad line region are visible, or Seyfert II galaxies, if the core is obscured 
by the torus [15]. Up to now, no Seyfert- like object with beamed emission has 
been observed. 



^ The spectral index a usually is determined between 2.7 GHz and 5 GHz. 



6 




Fig. 1. Scheme of an AGN with a black hole in the center and an accretion disk 
perpendicular to the direction of two jets along its rotation axis. The different incli- 
nation angles of the line of sight with respect to the jet for blazars, steep spectrum 
radio quasars and radio galaxies are indicated by arrows. Image adapted from [6]. 
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Fig. 2. AGN classification according to host galaxy, luminosity and inclination angle. 
2.1 GPS and CSS 



There is a substantial fraction of sources (~ 1/3 when selecting the strongest 
sources at 5 GHz) which shows, for frequencies above a certain turnover value, 
a steep radio spectrum from the compact radio core. If the turnover is in the 
MHz range, the source is classified as Compact Steep Spectrum source (CSS), 
if it is in the GHz range, the source is a CHz Peaked Source (GPS). These 
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sources arc significantly smaller than usual AGN: the GPS show linear sizes 
below 1 kpc, the CSS in the range of 1-15 kpc. GPS and CSS are not distinct 
source classes, since there is a continuous transition from the compact GPS 
to the slightly larger CSS. As pointed out in [9], the turnover frequency 
decreases with growing source size I as 

l^m OC I-'-'' . (1) 

The extremely high power of the GPS/CSS suggests a central engine similar 
to those of other AGN to provide sufficient energy. 

The source compactness can be explained by the assumption that jets get 
stopped by interactions with dense matter. The most common interpretation 
of the turnover observed for this source class is the model of synchrotron self- 
absorption, although free-free absorption on thermal electrons is not excluded. 
The sharpness of the maxima found imphes that the zone of the bulk radio 
emission is very small [16], presumably located at the outer end of the jet. 
GPS/CSS could be young states of radio-loud AGN evolving into larger radio 
sources. The population statistics of GPS/CSS and extended radio-loud AGN 
are in agreement with this picture, provided they stay in the compact state 
for a significant part of their lifetime. 

For neutrino production, this class is of particular interest if proton acceler- 
ation takes place within the central part of the jet. As a matter of fact, the 
dense matter surrounding the source provides an ideal target for pion pro- 
duction. Since jets get stopped in the interaction zone, the predicted neutrino 
production is isotropic or only slightly beamed. Thus, the inclination angle is 
considered to have only small effects on the expected neutrino flux on Earth 
from this sources. 

2.2 Radio-loud AGN with 100 kpc scale jets 

Two different jet morphologies have been observed for radio-loud AGN with 
extended jets correlated with the radio luminosity at 178 MHz [17]. These 
AGN can be distinguished according to their luminosity: the critical value 
of the luminosity is Liyg = 2.5 ■ 10^^ W/Hz (see Fig. 2), corresponding to a 
bolometric luminosity of 10"^^ erg/s [18,19]. High luminosity AGN are charac- 
terized by powerful jets extending far outside the host galaxy. The increasing 
jet luminosity at the outer end produces extended radio lobes and the so-called 
hot spots. These objects appear as Flat Spectrum Radio Quasars (FSRQ), as 
Steep Spectrum Radio Quasars (SSRQ) or as FR-II radio galaxies [17]. Low 
luminosity AGN, divided in BL Lacs or as FR-I radio galaxies, have fainter 
jets. They show decreasing radio emission with growing distance from the 
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central black hole and have no hot spots. 

This classification of radio-loud AGN assumes a typical value of the Lorentz 
factor 7 10 for the bulk motion in the jet [6] . This value can also explain 
the apparent superluminal motion of radio knots in blazar jets [20]. 

Radio-loud AGN show different appearances for different inclination angles. 
This is due to the relativistic Doppler boost of the emission from the jet and 
the obscuration of the inner core by the torus. For large viewing angles to the 
jet axis (close to 90°), the torus obscures the inner part of the AGN. Hence, 
the broad line region and the thermal continuum radiation of the accretion 
disk cannot be seen. In this case, the AGN is called a radio galaxy. 

If the opening angle of the torus is large enough, there is a range for the 
inclination angle, where the inner core is visible and the relativistic Doppler 
factor is smaller than 1. In this case the core is seen with its blue bump from 
the inner accretion disk and broad emission lines are present in the spectrum. 
The radio spectrum of these objects is still steep and lobe-dominated. Only 
high luminosity objects show this morphology, then appearing as SSRQ with 
a bright optical core and strong broad emission lines. The lack of similar 
low luminosity objects is still a matter of debate. Possibly, the inner core is 
obscured by the torus until the inclination angle is so small that the emission 
becomes Doppler boosted [19]. 

For very small inclination angles (< 12°), the jet radiation is Doppler boosted 
due to the relativistic motion of the bulk outflow towards the observer. These 
objects are characterized by a flat radio spectrum, strong variability and po- 
larization. FSRQ and BL Lac objects can be combined into the blazar class 
since both are characterized by strong beaming effects. The flat radio spectra 
of blazars can be explained by the dominance of the boosted flat-spectrum 
core over the non-boosted steep-spectrum radio lobes and will be used as se- 
lection criterion for blazars. A flat spectral index in radio is an indication of 
optical thickness [13]. 

BL Lac objects are low-luminosity objects with a FR-I type jet, while FSRQ 
are considered as the high luminosity objects (FR-II jet). The BL Lac objects 
are commonly divided into high-energy cutoff BL Lac (HBL) and low-energy 
cutoff BL Lac (LBL), referring to the maximum energy of the electron syn- 
chrotron spectrum. HBL are relatively weak in radio flux, strong in X-ray 
flux and bolometrically less luminous than LBL. A distinction between these 
classes is usually made by taking into account whether the radio-to-X-ray spec- 
tral index (with F(u) oc ) is bigger (LBL) or smaller (HBL) than -0.75 
[6]. The high energy 7 emission of HBL and LBL differs: all confirmed AGN 
TeV sources are HBL^ [21]. In contrast, the stronger AGN GeV (EGRET) 

^ However, there are indications of TeV radiation from BL Lacertae (LBL) and 
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sources consist of LBL and FSRQ, and only 2 of the 6 confirmed TeV 7 AGN 
have been detected by EGRET, with only moderate fluxes [22]. 



3 Selection of neutrino source candidates 

Accelerated proton interactions with ambient photons or matter lead to neu- 
trino production through the reactions: 

— > nn ' 
pp^TT^ + N 

and similar reactions for neutrons. These processes are always accompanied 
by neutral pion production. Neutrino and gamma ray fluxes are expected to 
be of the same order of magnitude. Since in most scenarios, proton and tar- 
get photon spectra fall steeply with energy, higher resonances and multi-pion 
production represent a small correction to pion production. The acceleration 
of the protons is assumed to be due to shock acceleration which may take 
place in the relativistic jet [20,23] or in the accretion disk [24,25]. Due to low 
plasma density in AGN jets, the pj interactions are likely to be dominant over 
pp [20]. This is also indicated by the lack of absorption lines in AGN X-ray 
spectra [26]. An exception from this consideration is given, if pp interactions 
take place at the inner edge of the accretion disk [25]. 

If acceleration takes place in the disk, electromagnetic cascades initiated by 
photons from neutral pion decay would lead to a non-thermal X-ray spectrum. 
However, measured AGN X-ray spectra and the diffuse X-ray flux show the 
dominance of a thermal flux peaked around 100 keV [27]. This discrepancy 
might be solved by assuming that the non-thermal X-ray spectrum contributes 
only 30% to the measured X-ray flux [26]. That implies that hadronic inter- 
actions are not dominant. Hence, the resulting neutrino spectrum would be 
reduced by the same factor. Due to this discrepancy, the disk model is disfa- 
vored as a dominant process for neutrino production. 

For the jet as well as for the disk, the theory of diffusive shock acceleration 
suggests a power law spectrum with a differential spectral index a ^ 2 for 
protons [28,29,30]. Neutrinos and photons are produced with the same index 
if multi-pion production is not dominant. However, the energy of photons 
produced by neutral pion decay may be redistributed to photons of lower 
energy by synchrotron pair cascades. The photon escape energy depends on 

M87 (misaligned BL Lac or FR-I) weakening the association. 
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the optical depth. Hence, the spectra of neutrinos and photons can differ 
considerably, while the bolometric luminosity remains correlated [31]. 

Hadronic jet models have been proposed to explain the high energy 7 emission 
from the jet. A comparison to purely leptonic models can be found in [20]. 
Purely leptonic models usually explain the high energetic (TeV) radiation by 
inverse Compton scattering of photons. Models in which the photons originate 
as synchrotron radiation are classified as Synchrotron Self Compton (SSC); 
those in which the photons are ambient are classified as External Compton 
(EC) . The EC model requires high densities of external photons to be Compton 
up-scattered to higher energies. These photon densities then constrain the 
escape energy of photons to values below 300 GeV. Thus, the leptonic EC 
model cannot explain the observed TeV sources [20]. In the SSC models, the 
ratio between GeV and radio luminosity is constrained, predicting a larger 
number of GeV blazars than observed by EGRET. Additionally, keV and 
TeV fluxes are correlated. Many observations confirm these correlations and 
apparently favor the SSC model. However, recent observations indicate the 
existence of orphan flares for lES 1950+650 [32] and for Mrk 421 [10,33], 
where the TeV flux flares but the X-ray flux does not. 

In hadronic models, the emission of high energetic photons in the jet is not 
assumed to be directly from photohadronic processes. A non-negligible opti- 
cal depth makes TeV photons cascade to lower energies. Thus, an additional 
mechanism for high energy gamma production is necessary. In proton initiated 
cascade models (PIC), it is assumed that electrons scatter on some target pho- 
tons acquiring higher energies. In the SS-PIC (self synchrotron PIC) model, 
the target photons are those emitted in the cascades. The EC-PIC (external 
Compton PIC) assumes external photons to be scattered up to higher energies 
by the inverse-Conipton process [20]. 

Another possibility to explain observed high energy 7 emissions is given by the 
Synchrotron-Proton Blazar model proposed by Miicke and Protheroe [34]. It 
assumes relativistic protons to emit synchrotron radiation at higher energies 
than the one from gyrating electrons. This model explains the observed double 
hump energy spectrum and predicts neutrino emission spectra for HBL and 
LBL [35] . The emission of TeV 7-rays via synchrotron radiation requires pro- 
tons to reach extremely high energies, which is possible only if the interaction 
rate is sufficiently low. Hence, in TeV emitters and pp interactions must 
be less frequent than in blazars without intense TeV photon emission. This 
results in predictions of neutrino fluxes for HBL being 6 orders of magnitude 
less than for LBL. 

In summary, to select source candidates, we assumed the common origin of 
neutrinos and photons from pion decay and the classification described before. 
Since the optical depths of AGN are unknown, escaping photons can be of 
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significantly lower energy. This is taken into account by the selection of source 
candidates at various energies where data are available (radio, IR, optical, keV, 
GeV, TeV). Moreover, we assume that the observed photon flux is proportional 
to the TeV neutrino signal. While this selection does not depend on detailed 
AGN models, other explicitely model dependent selections are possible, e.g. 
based on the concept of the jet-disk symbiosis [13,18] used in [36]. 



4 Selection of neutrino candidate sources from catalogs of AGN 
classes 

We consider catalogs of AGN classes and apply some selection rules, described 
in Tab. 1, in order to obtain 'statistically complete catalogs' of sources in a 
well-defined part of the sky above a fiux threshold. In this selection we did 
not consider the variability in time of emissions from AGN. As a matter of 
fact, most of the available photon data were not collected at the same time as 
the considered AMANDA data set. Future multi- wavelength campaigns and 
simultaneous collection of photon and neutrino data will allow approaches 
that also account for time dependencies. 

To all the considered catalogs we applied the requirement that the galactic 
latitude b is larger than 10° (except for radio-weak quasars, where the catalog 
itself requires |6| > 30°). This cut excludes the galactic plane in order to avoid 
biasing AGN samples with galactic sources. Moreover, a minimum declination 
of 10° is required since we look for neutrino induced upgoing muons. We also 
consider a threshold on the distance corrected fiux, to exclude intrinsically 
weak nearby sources which might otherwise migrate into the catalog due to 
their distance. The formula used for the correction is given in appendix A. 
We stress that some blazar sources have been removed from our catalogs by 
a cut on the distance corrected flux. The sources that were removed only due 
this cut are listed in Table B.ll. Among these there are some well-known 
nearby AGN, e.g. NGC 1275 would be in the IR-blazar and in both keV- 
blazar samples and Mrk 501 would be in both keV samples. The occurrence of 
these sources in various samples confirms our assumption that these sources 
migrate into the samples due to their proximity and should not be considered 
as generic sources. These nearby sources with such high fluxes would bias a 
stacking analysis, hence they should be analyzed as individual sources [1]. 

In the following sections we describe with more details the catalogs for the 
considered classes of sources in Tab.l to which we apply selection rules. Then 
we consider for which of these selected sources there exist measurements of 
the photon flux that we assume to be correlated to neutrino fluxes and give 
our lists of candidate neutrino sources for the various classes. 
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Sample 



Flux /luminosity cuts 



Coordinate cuts 



Further cuts 



Blazars 
CSS 
GPS 
FR-I 

FR-II 

QSO 



F^GHz > 0.8Jy, Fl^fjl > 1.5 Jy 
Fi78MHz > lOJy, logPnsMHz > 26.5 
F5GHZ > IJy, iogP^GHz > 26.5 
FnsMhz > lOJy, 
Ll78MHz < 2.5 • 1025 W/Hz 

Fi78Mhz > lOJy, 
Ll78MHz > 2.5 • 1025 W/Hz 
B < 16.16, U-B < -0.44 



> 


10°, 


\b\ 


> 


10° 


> 


10°, 


\b\ 


> 


10° 


> 


10°, 


\b\ 


> 


10° 


> 


10°, 


\b\ 


> 


10° 



a > -0.5 with F^, oc 1^°= 

Lin. size< 20kpc 
0.4GHz< Fmax < 6GHz 



d > 10°, |6| > 10° 



S > 10°,|6| > 30° 



Table 1 

The selection criteria defining catalogs of the different source classes. Fi, stands 
for the photon flux at a certain frequency u in Jy, F^=°^ is the distance corrected 
flux to a redshift of z = 0.1, logP is the logarithmic power in W/Hz and L is the 
luminosity. The declination is labeled as 5 while b stands for the galactic latitude. 
For the blazars, the spectral index a is determined between 2.7 GHz and 5 GHz. 
For CSS selected at 2.7 GHz, a back-extrapolation of a steep spectrum component 
to 178 MHz was used instead of the measured 178 MHz flux. U is the magnitude 
of the UV flux (A = 300 nm-400 nm), B the magnitude of the flux at A = 400 
nm— 550 nm (blue). 



4-1 Blazars 



As explained in Section 3, blazars are characterized by a flat radio spectrum, 
which we use as the selection criterion here to create a 'blazar hst'. For radio 
sources with high fluxes at 5 GHz (A = 6 cm), Kiihr et al. [37] worked out a 
complete catalog of radio sources with flux larger than about 1 Jy ^ . A second 
version also covers sources below this threshold, using data available up to 
1981 [38]. In the Northern sky, this catalog contains the "strong surveys" Sl- 
S5 with a threshold of 0.8 Jy and below. The selection rules select from the 
blazar catalog about 150 sources. 

Since blazar jets are in the direction of the observer, blazars are interesting 
candidate sources of neutrinos that could be produced in the jet. Correspond- 
ing to different assumptions on the optical depth of blazar jets, the neutrino 
flux is assumed to be proportional to the photon flux that are measured by 
various experiments in different energy ranges: HEGRA and Whipple in the 

^ IJy = 10~26ws-i]:n-2^ older papers also called f.u. for flux unit. 
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TeV one, EGRET in GeV, HEAO-A and ROSAT in keV or IRAS in infrared. 
The ordering of the sources is done according to the flux vahie in the relevant 
energy range. For the resulting samples, the relative source strengths at the 
selection energies are displayed in Fig. 3(a). The flux of the strongest source 
is normalized to 1. 

For the sources in our list we look for catalogs of the measured photon flux in 
various bands to order them according to their photon emissions. 



IR-loud blazars: the Infrared Astronomical Satellite (IRAS) provides a sur- 
vey at wavelengths 12/xm, 25//m, 60//m and 100//m. The Faint Source catalog 
IRAS F, covering most parts of the sky outside the galactic plane, includes 
173,044 sources [39]. We find in total 12 sources in the IRAS F catalog that 
satisfy our selection criteria and match sources in our original 'blazar list'. 
The IR selected sample is given in Table B.l. 



X-ray blazars: in the X-ray regime, there are only a few all-sky surveys. 
We make the hypothesis that soft and hard X-ray emissions are positively cor- 
related with TcV neutrino production. For hard X-rays (0.25-25 keV) there 
exist only the HEAO-A measurements performed in the 1970s with low sensi- 
tivity and relatively poor angular resolution. This resulted in the IH catalog 
[40] containing 842 sources. The X-ray satellite ROSAT performed an all-sky 
survey in the 0.2-2 keV band, with high sensitivity and good angular resolu- 
tion in the range of a few tens of arc-seconds. The strongest sources found in 
the all-sky survey are listed in the IRXS catalog [41], which contains in to- 
tal 18,811 sources. For these antipodal satellites, two separate candidate lists 
are constructed from IRXS and IH sources that match our 'blazar list'. This 
resulting selections are hsted in Table B.3 for HEAO-A and in Table B.2 for 
ROSAT. 



GeV blazars: in the GeV range, the most sensitive data are those col- 
lected by the EGRET experiment on board the CGRO satellite [42]. The 
third EGRET catalog (3EG) [22] contains in total 271 sources of which 66 
have been identified as blazars and additional 27 sources with lower confi- 
dence. We select EGRET sources that are rehably identified with a blazar in 
our list (class 'A' in the 3EG catalog), while weakly identified sources (class 
'a' in the 3EG catalog) are accepted if an independent analysis confirms this 
identification [43]. The sources are sorted according to the maximum flux de- 
tected by EGRET. This yields the source subsample observable by AMANDA 
presented in Table B.4. 
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Wc include a further sample (see in Tabic B.5) selected among the unidentified 
EGRET sources (these have not yet been identified with sources detected at 
other wavelengths). Although a large fraction of these sources is likely to 
be galactic, some may be extragalactic. To select presumably extragalactic 
sources, aside from the usual galactic plane exclusion, we check that sources 
do not show features typical for galactic ^ sources. Selected sources are listed 
in Table B.5. These sources show a high variability, even higher than that of 
the identified AGN, and a large energy output that can be explained, even if 
they are also AGN. 



TeV blazars: up to now, only 6 blazars have been firmly detected emitting 
TeV photons, of which 5 are in the field of view of AMANDA [21]. All TeV 
blazars show extreme variability in flux, such that a flux ranking becomes 
impossible. The sources all have a quite weak radio flux and do not fulfill 
our radio selection criteria for blazars. Furthermore, due to absorption by the 
IR background, only nearby sources can contribute to the TeV photon fiux. 
These conditions prevent the use of our selection scheme for TeV blazars. But 
since in optically thin sources. TeV photons are strictly correlated to TeV 
neutrinos, we consider all observable TeV blazars, listed in Table B.6, for the 
source stacking analysis. 
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Fig. 3. Relative strength at the selection energy of strongest sources in blazar sam- 
ples (a) and in non-blazar samples (b). 



^ Three different types of presumably galactic sources have been considered: two 
are concentrated on the galactic plane with |6| < 5°, while the third is associated 
with the Gould belt [44]. The latter form the so-called Local Gamma-Ray Popula- 
tion characterized by constant fluxes below 2.4 • 10~^s~^cm~^ and a soft spectrum 
(photon spectral index 7 2.45); they reach galactic latitudes up to 30°. 
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12 GPS and CSS 



For defining a catalog of compact sources, we follow the selection of O'Dea 
[9]. The catalog contains former selections of CSS by Fanti et al. [45] and 
of GPS by Stanghcllini et al. [46]. Since CSS and GPS are characterized by 
a maximum in the radio flux at a size-dependent frequency, the selection of 
these sources is done at several frequencies: 178 MHz, 2.7 GHz and 5 GHz 
depending on the peak location. 

From the original CSS samples, all sources have been removed which do not 
fulfill all the selection criteria in the original paper [45]. The same has been 
done for CSS and GPS sources when the requirements of the more recent 
analysis [9] are not fulfilled. Additionally, we put a cut on the absolute power 
at 5 GHz, similar to the cut on the power at 178 MHz for the CSS sources 
indicated in Tab. 1. The typical spectra of sources arc plotted in Fig. 4. In 
total 40 sources are selected of which 9 are GPS, 23 are CSS and 8 sources 
fulfill selection rules for both classes. 




178 MHz 5 GHz ^ 178 MHz 5 GHz ^ 178 MHz 5 GHz ^ 



Fig. 4. Schematic plot of radio spectra of sources fulfilling the CSS and GPS selection 
rules. On the left: a typical spectrum of a CSS selected at 178 MHz. In the middle: 
a CSS spectrum with the maximum near 2.7 GHz. On the right: a GPS spectrum 
selected at 5 GHz. For the sources selected at 2.7 GHz, the back-extrapolation of 
the steep spectrum component to 178 MHz is relevant for the selection. 

Since CSS and GPS were selected at three different radio frequencies, we 
cannot use the intensity in the radio. Since there are only few X-ray or IR 
data for these sources, the only possibility to get a suitable source ranking is 
to use data in the optical band taken from [9] . The underlying hypothesis we 
make is that the TeV neutrino fiux is proportional to the optical fiux. Our 
resulting neutrino source catalog is displayed in Table B.7. 

4.3 FR-I and FR-II-radio galaxies 

As the distinction between FR type I and type II radio galaxies is made at 
178 MHz, it is adequate to select radio galaxies from observations at this 
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frequency. Data at this frequency can be found in the 3C catalog in its revised 
version (3CRR) [47]. The classification into FR-I and FR-11 galaxies depends 
on the availability of the redshift for the source and is taken from [48]. The 
selection results in a catalog of 24 FR-I and 122 FR-II galaxies. Additionally, 
16 compact sources fulfill the selection. They belong to the CSS source class 
as discussed in 4.2 above and are removed from the FR-I and FR-II samples. 



For the selection of the neutrino candidate list, the radio fiux at 178 MHz 

is assumed to be proportional to the TeV neutrino fiux. Hence, we use these 
fiuxes to sort the samples of radio galaxies. The results are listed in Table B.8 
for FR-1 and in Table B.9 for FR-II galaxies. Since FR-I/II radio galaxy jets 
do not point towards us, the observable neutrino fiux produced in the jets is 
expected to be weaker for these sources than for blazars. Nonetheless, neutrino 
production in the disk is expected to be isotropic compared to the case of jet 
emission. This enhances the observation probability since it is not required 
that the observer is in the direction of the jet. 



4-4 Radio-weak Quasars 



A complete catalog of quasars selected at optical to UV wavelengths is given 
by the Bright Quasar Survey (BQS) [49]. A quasar is identified by an excess 
of UV fiux. The catalogs presented by Sanders et al. [50] contain in total 109 
sources, of which 59 sources are radio-weak quasars in the sky viewable by 
AMANDA. 



Mult i- wavelength investigations of the BQS sample show a second maximum 
in the photon spectrum at IR wavelengths for some of the quasars [50]. Hence, 
in the usual assumption that photon fiux features are correlated to neutrino 
production at the source, we include in our list the radio- weak quasars selected 
according to their photon fiux at 60 fim. The full quasar catalog selected from 
the BQS according to the IR fiux is given in Table B.IO. 



For the resulting samples of CSS/GPS, FR-I/-II radio galaxies and radio-weak 
quasars, the relative source strengths at the selection energies are displayed 
in Fig. 3(b). The fiux of the strongest source is normahzed to 1. 
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5 Optimization of the number of sources and the bin size 

5. 1 The statistical optimization procedure 

The different hypotheses on the neutrino flux from the sources in our catalogs 
determine the relative source strength within each class but leave the normal- 
ization free. We have normalized the signal by relating the sensitivity of the 
point source analysis to the signal of the strongest source in each sample. We 
consider here data taken by AMANDA during the year 2000 [1] . 

In the point source analysis, a search was performed for a statistically signif- 
icant excess from any direction. For this purpose, any position on a grid of 
rectangular search bins is evaluated by considering the background expecta- 
tion calculated from the average number of events in the corresponding zenith 
range. No point sources were detected; the observed excesses were equivalent 
to those in a randomized distribution. In the considered data set, we find by 
multiplication of the sensitivity^ with the effective area, the lifetime of the 
analysis [1] and the bin efficiency for a circular search bin of 3.5° radius that 
the sensitivity (90% C.L.) corresponds roughly to three neutrinos in a 3.5° 
bin. This motivates the choice of the normalization of the neutrino flux for 
the most intense source that would produce a neutrino signal Si <3 events in 
the considered AMANDA data set. Hence, for this source stacking analysis, 
we optimized the number of sources using 5*1 = 1,2,3 events (also using the 
bin efficiency for 3.5°). Significantly higher signal normalizations would imply 
that the most intense source must be correlated with the position of a few 
spots, where an excess has been observed. These excesses are most probably 
due to statistic fluctuations. Since we do our selection independent of the re- 
sults from that data set, such a coincidence is not assumed. On the other hand, 
for 5*1 < 1, the expected signal is too faint to result in a notable contribution 
in the source stacking analysis. For the small numbers considered, Poisson 
statistics has to be applied. With our normalizations, the assumed neutrino 
fluxes for the sources in our catalogs agree with the limits of the current diffuse 
analyses [51,52]. 

Since the sensitivity of the AMANDA-II detector is almost independent of the 
declination, in the optimization procedure we approximate the background to 
be constant in the considered declination range. A total amount of 699 events 
in the year 2000 neutrino sample results in a background of 111.25/sr for that 
year. 

A given number of sources Ngrc corresponds to a mean number of expected 
signal events, S, and a mean number of expected background events, BG. 

^ We define the sensitivity as the average upper limit in case of no signal. 
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Using Poisson statistics, we optimize for the median significance which will be 
reached in 50% of the experiments being performed under these assumptions. 
First, we calculate the median number of events n under signal hypothesis 
as a function of the number of sources to be included. Then we calculate the 
Poisson probability to observe at least n events with the assumption of pure 
background. Rescaled to the corresponding number of standard deviations of 
a Gaussian, 5, we get the median significance as a function of the number of 
sources to be included. The optimum number of sources, as determined by 
this median significance procedure, corresponds to the highest value of 5. 

The same procedure is apphed to the size of circular search bins, varying the 
signal according to the point spread function (PSF) and the background ac- 
cording to the area covered by the search bins. The PSF is evaluated for the 
most hkely case of an E~'^ spectrum. Since the PSF is not completely indepen- 
dent of the dechnation, we use the PSF averaged over all dechnations. In Fig. 5 
the average PSF obtained by considering an isotropic source distribution over 
the sky with 5 > 10° is plotted for various spectral indices. The declination 
dependence of the median of the PSF is shown in Fig. 6 for an spectrum. 
The resolution slightly improved for larger values of the dechnation. 
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Fig. 5. Cumulative Point Spread Function after quality cuts for various spectral 
indices in year 2000 sample. 



The optimum number of sources and the optimum bin size are not independent 
of each other. Thus, we follow an iterative procedure. First we evaluate the 
number of sources with rough estimates for signal and background. Then we 
optimize the bin size for the given number of sources, and finally, we check 
the number of sources again. 

The results for the optimum values have been checked by a second procedure 
which will be called maximum observation probability procedure. A minimum 
significance 6 was predefined and the probability to observe at least this sig- 
nificance is optimized. As the expected signal is rather small, a value of 3a 
was chosen for 5. 



19 



a 2 



0.5 




AMANDA (2000) 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Sin(DEC) 



Fig. 6. Median PSF of AMANDA for data collected in 2000 as a function of 
sin{DEC). 
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Fig. 7. Median significance as expected under the considered hypothesis as a func- 
tion of the number of GeV blazars to be stacked. The highest value of the median 
significance corresponds to the optimum number of sources. 

5.2 The resulting samples 



We use the total number of events and the PSF from the year 2000 point 
source sample [1] for the optimization to determine samples from our catalogs 
to be analyzed for a cumulative neutrino fiux. Our procedures may be applied 
to other data sets easily. The resulting parameters for year 2000 data are listed 
in Table 2. 



In Fig. 7, we plot the median significance as a function of the number of 
EGRET blazars to be included for signal normalizations 5*1 ranging from 1 to 
4. The highest normahzation of 5"! = 4 signal events from the strongest source 
is unlikely, as it implies a median significance above 3a for the single source 
search. Hence, Si — A has not been taken into account when evaluating the 
optimum. 

For most source classes, we find an optimum of about 10 sources. With one 
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exception, the source stacking analysis is more sensitive in testing the corre- 
sponding hypotheses than the point source analysis of single sources. Usually, 
the observed peaks are asymmetric: if more sources than the optimum num- 
ber are included, the expected significance falls only slowly. On contrast the 
significance falls quite steeply if fewer sources are included. 

The exceptions arc found for FR-1 and for FR-II radio galaxies, where the 
results are completely different. As there is a luminous FR-I galaxy in our 
neighborhood, M87, the whole FR-I fiux would be strongly dominated by M87. 
The detection of TeV photons from this non-blazar AGN by HEGRA [53] and 
H.E.S.S. [54], theoretically predicted in [55], also supports the uniqueness of 
this source. Thus, the optimization suggests just to analyze the neutrino fiux 
from that source only. In other words, the analysis of the strongest source only 
is more sensitive to our hypothesis than source stacking for the corresponding 
catalog. However, for two reasons we decided to analyze additionally a FR-I 
sample without M87. First, there is always the possibility that a single source 
does not contribute to the neutrino flux due to a reason not considered at 
selection. Second, in the opposite case, if M87 would be identified as a neutrino 
point source, the cumulative signal of the other sources contains information 
whether the M87 neutrino fiux would be a specific feature of that source or 
a general characteristic of FR-I galaxies. Thus, we analyze a FR-I sample 
without M87 additionally. Then we find an optimum value of 20 sources to be 
included. 

The situation of FR-II galaxies is opposite to that of FR-I galaxies. There are 
a lot of sources at higher redshifts with similar luminosities. This results in an 
monotonically increasing sensitivity while including more sources. Eventually, 
the coverage of the sky by search bins reaches a significant fraction of the sky 
and the techniques of point source analysis become less suitable to analyze 
the signal. If the energy spectrum of the signal differs from the atmospheric 
background, diffuse techniques searching for global excess of high energy neu- 
trinos should be more sensitive to FR-11 galaxies. However, since point source 
analysis and diffuse analyses are using complementary techniques and a point 
source analysis can give more information about the origin of the excess, a 
source stacking analysis of FR-II galaxies may also be useful. Since our op- 
timization fails, we have to define the most suitable value for the number of 
FR-11 sources to analyze. For that, we use a local saddle point at 17 sources 
for the cutoff where the sensitivity grows only marginally if further sources 
are included. 

The unidentified EGRET blazars show nearly constant values for observation 
probability and median significance between 15 and 35 sources. 

The hard X-ray sources are restricted to three sources due to the poor sensi- 
tivity of HEAO-A, an experiment of the 1970s. 



21 



We cross-checked the procedure based on the mean significance to that based 
on the maximum observation probabihty. The deviations between the resuhs 
from the two procedures are within the corresponding uncertainties of the 
determination of the optimum. 



For the optimum bin size, we find resuhs between 2.6° for the unidentified 
GeV sources and 3.0° for the single source M87. As expected [56], the bin size 
decreases with the size of the sample. 



source class 


ATMS 

src 


NDP 

src 


listed in 


Bin size [deg.] 


IR blazars (IRAS) 
keV blazars (ROSAT) 
keV blazars (HEAO-A) 

GcV blazars 
unidentified GeV sources 
TeV blazars 
GPS and CSS 


lit? 

8±1 

3±0 
8± 1 
22 ±5 

5* 
8± 1 


111? 

9± 1 

3±0 
9± 1 
20 ± 5 

5* 
8±2 


table B.l 
table B.2 
table B.3 
table B.4 
tabic B.5 
table B.6 
table B.7 


2.8 ±0.2 

2.8 ±0.2 

2.9 ±0.2 
2.8 ±0.2 
2.6 ±0.2 
2.8 ±0.2 
2.8 ±0.2 


FR-I radio galaxies 
FR-I without M87 
FR-II radio galaxies 
Radio-weak quasars 


1±0 
20^3 

122 


1±0 
19 ±4 


table B.8 
table B.8 
table B.9 
table B.IO 


3.0 ±0.2 
2.6 ±0.2 
2.6 ±0.2*** 
2.8 ± 0.2 



Tabic 2 

Resulting parameters for the source stacking analysis. 

and stand for the optimum number of sources determined with the me- 
dian significance procedure and, respectively, the maximum observation probability 
procedure. 

* All sources included without optimization, see above . 

** This value was taken for the analysis. Only a saddle point. The observation prob- 
ability increases again, when including all sources (— diffuse analysis). See Sec. 5.2. 
*** Bin size evaluated for 17 sources. 

5.3 The relation between flux distribution and the optimum number of sources 

As described above, the optimization of the number of sources can lead to 
two special cases, where either the source stacking analysis degenerates to 
a point source analysis of a single candidate source, or to a diffuse analysis 
where the arrival direction of the neutrino becomes irrelevant. If the flux falls 
only slowly from the strongest source to the weakest, as is the case for FR-II 
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radio galaxies, then we cannot find an optimum number of sources. So the 
probabihty to observe a signal always increases when more sources are added. 
In the other extreme, as was found for FR-I radio galaxies, if the strongest 
source is expected to contribute most of the flux and the other sources are 
substantially weaker, then the sample of sources to be stacked degenerates to 
a single source. Here we investigate the conditions on the flux distribution for 
these degeneracies to occur. We consider the flux F of the A^*'^ source, F{N). 
Motivated by the nearly linear decrease of F{N) in the double logarithmic 
plots in Fig. 3, we assume a power law, 

F{N) oc AT" . (2) 

For a > —1, the finiteness of the total flux requires a cut-off or a steepening. 

For different values of a, we evaluate the optimum number of sources and 
determine the range of values for a, where the two degeneracies occur. The 
maximum number of considered sources is set to 100. This corresponds to 
a coverage of about 10% of the sky viewable by AMANDA, assuming a bin 
size of 2.6° . For significantly more sources, the sky get densely covered by 
candidate sources and the directional information cannot be used and diffuse 
methods are more useful. 

For a slowly falling flux distribution with a > —0.65, we cannot flnd an 
optimum number and all sources have to be included. In the range —0.65 > 
q; > — 2, a non-degenerate optimum number of sources is found and for a < 
—2, only the strongest source has to be selected. The thresholds depend on 
the detector parameters, although for a > —0.5, S/\/BG diverges, i.e. the 
degeneracy occurs generally for these values. An improved angular resolution 
can explore the range —0.65 < a < —0.5 for a source stacking analysis. 

A fit of a on the flux distributions of our samples results in o; = —0.6 for 
FR-II radio galaxies and for IRAS blazars, a 2-point index a = —4.1 for the 
two strongest FR-I radio galaxies and —0.7 > a > —1.1 for the other samples. 
For all samples, the results of the optimization are in agreement with our 
classification according to a. 



6 The results for the year 2000 

The sources we selected are hsted in Table B.4-B.10 in the appendix. We ap- 
ply the source stacking method to the data collected by AMANDA in the 
year 2000 using the event sample from the previous point source analysis [1] 
corresponding to a detector lifetime of 197 days. The selection of well recon- 
structed events was optimized for high sensitivity to point sources. In Fig. 8, 
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the reconstructed arrival directions of the individual events of this data set 
are shown in a sky map. The zenith distribution of these events is displayed in 
Fig. 9 in comparison with the simulated angular distribution of atmospheric 
neutrinos. Since the event selection was not optimized for purity, the excess of 
measured events close to the horizon is expected as the contribution of poorly 
reconstructed atmospheric muons. The detector response with this event se- 
lection was simulated for various spectral indices. The energy spectrum folded 
with the Earth absorption factor and with the detector response is shown in 
Fig. 10. A uniform source density on the Northern Sky was used in the simu- 
lation, corresponding to a mean zenith angle of 33.4°. For an power law 
spectrum, as expected from Fermi acceleration, that means 90% of the signal 
events are between the AMANDA energy threshold at 50 GeV and 350 TeV. 




-90° 



Fig. 8. Reconstructed arrival direction of 1557 selected events in the year 2000 
sample [1]. The events below 5 = 0° are dominated by atmospheric muons. 
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Fig. 9. Distribution of the zenith angle in the considered data set in comparison 
with the expectation of atmospheric neutrinos from simulations. 

Results obtained on a more complete set of data collected in the years 2000- 
2003 will be reported in [3] . Moreover, this method and the list of sources we 
derived will serve as the starting point for stacking analyses performed with 
the IceCube neutrino telescope [57] which is currently under construction. 
AMANDA will be integrated into IceCube. The source stacking method will be 
adapted to the better performance of the larger detector in terms of effective 
area and angular resolution. Application of this method for other neutrino 
telescopes [58,59,60,61,62] is conceivable. 
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Fig. 10. Energy spectra folded with the response of the AMANDA detector for 
hypothetical source spectra with the spectral indices —1, —2 and —3. 

In this analysis, only the cumulative result for a class of sources is evaluated, 
not the signal for individual sources. The AMANDA collaboration follows a 
strict blindness policy, i.e. an analysis has to be developed blindly with respect 
to the data. This prevents statistical fluctuations to affect the final steps of 
the analysis. In the context of the source stacking analysis, the selection of 
sources was done independently of the AMANDA data itself. No signal was 
evaluated until the source samples and all analysis parameters were fixed. 

The number of observed events in our search bins has to be compared with the 
background, mainly atmospheric neutrinos. Since the background is indepen- 
dent of the right ascension, for each position in the sky it can be evaluated from 
the events in the same declination band. If two sources of the same sample are 
very close to each other, their search bins may overlap. Then, the background 
expectation is corrected for the overlap and events in the overlapping area 
contribute only once to the cumulative signal. 

The probability to observe the measured number of events under the hypothe- 
sis of pure background is evaluated by Poisson statistics. The correctness of the 
analysis was tested with two procedures. First, to test the correct evaluation 
of the significance, a collection of data sets with randomized right ascension 
was created and the significance of the signal from the stacked sources in 
these randomized data sets was determined. As expected for randomized data 
sets, the significance follows a Gaussian with a width of one centered at zero. 
Additionally, hypothetical source lists with random source positions were eval- 
uated using the original data set. The results of the point source analysis of 
the considered data set [1] and the unfolded neutrino spectrum obtained from 
that data set [8] suggest that the data sample is strongly dominated by at- 
mospheric neutrinos. If this holds, a difference of the significance distribution 
to that of randomized data sets, which follows a Gaussian (see above), would 
indicate that our assumption of a fiat background in right ascension is not 
correct. Also in this case, the observed significance distribution follows our 
expectation. The significance distribution for both tests is shown in Fig. 11. 
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Fig. 11. Significance distribution for selected sources in scrambled data sets (a) 
and for samples of 10 sources with random source positions evaluated with the 
unscrambled data set (b). 

No statistically significant signal was found, in this limited data sample, above 
pure background; only upper limits could be derived. To calculate these hmits 
we used Feldman-Cousins [63] confidence intervals and Poisson statistics. Sys- 
tematic errors, estimated to be at the ±30% level, were not taken into account. 
The measured event rates, the background expectation and the resulting up- 
per limits in terms of event counts and integral neutrino fluxes above 10 GeV 
are presented in Table 3 and visualized in Figure 12. A dedicated analysis of 
systematics for a multi-year dataset is currently in preparation. For most of 
the declination range, the sensitivity (90% C.L.) of the point source analysis 
for the considered data set is about 2 • 10~^cm~^ s~^ [1]. A comparison with 
the hmit per source obtained from stacking (see Table 3) shows the progress 
reached by the source stacking method. 



7 Summciry 



We have performed a systematic classification of AGN. This is the basis of a 
source stacking analysis of TeV neutrinos with the AMANDA neutrino tele- 
scope. Neutrinos are assumed to be produced in optically thick jets or accretion 
disks. The optical depths of these sites and consequently the average energy 
of photons leaving the source are unknown. We therefore select the sources 
according to the photon fiux at different energies. Most of the available photon 
data were not collected at the same time as the considered AMANDA data 
set. Thus, the time variabihty of AGN was not considered in this analysis. 
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Table 3 

Results for the year 2000 data: Number of sources Nsrc, measured number of 
events N°^^ , the corresponding background N^^ and the 90% C.L. limits on the 
event counts (num) and on the integral flux for an E^'^ spectrum above 10 GeV 
iflim) in units of 10~^cm~^ s"^. fiim/Nsrc represents the limit per source. 
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Fig. 12. Limits on the neutrino flux from an average source in the generic AGN 
classes using data from the year 2000. The average sensitivity to single point sources 
is indicated by horizontal lines for this data set [1] a well as for a 4 year data set 
[3]. 

With the hypothesis that the TeV neutrino flux is proportional to the corre- 
sponding photon flux, we have optimized the number of sources to be included. 

For most source classes, an optimum of 8 — 12 sources to be stacked has been 
found. For these, the source stacking analysis is more sensitive than the point 
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source analysis of a single source. Exceptions have been found for FR-I and 
FR-II radio galaxies. While FR-I galaxies are dominated by the single source 
M87, there are too many FR-II galaxies with a similar flux in order to converge 
at a reasonably small number of sources to be stacked. 

For the data of the year 2000, no significant deviation from the background 
expectation was found. In a next step, the stacking method will be applied 
to several years of AMANDA data. Optimizing the method to an angular 
resolution of less than 1° and better background rejection, it will be adapted 
to the 1 km^ IceCube array [57]. Starting with the selection of southern sky 
sources, it is also applicable to other neutrino telescopes [58,59,60,61,62]. 
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A The distance correction of fluxes 



The distance dependence of a measured flux F,, at the observer's frequency u 
of a source with luminosity L,^/ is given by 



(A.l) 



with the luminosity distance di and the frequency at the source i/' — + z). 
The luminosity distance as a function of redshift is given by 



di{z) = < 



^^^sm{\n,f-'-I{z)} , ifOfe<0, 



(A.2) 



see [64]. I{z) is given by 



I{z) = J [{1 + z' f ■ (1 + Vt^z') - Vtxz'{2 + z')\ ' dz' , (A.3) 





-0.5 



with the normalized curvature flk- We used = 0.27, Q\ — 0.75 according 
to [65] and h — 0.71 according to [66]. 

For a power law spectrum with spectral index a, the redshift results in 

L,,du' = L,.(i+,) • (z + 1)"+Mi/ . (A.4) 



Finally, the correction to 2; = 0.1 for a source at redshift z is given by 



?z=0.1 



F- 



di{z) 



z+l 



di{z = 0.1)J VO.1 + 1 



a+l 



(A.5) 



B Source catalogs for AMANDA 
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Table B.l 

Strongest IR sources in blazar catalog. F^ghz stands for the radio flux at 5 GHz, 
while the mean IR flux at 12, 25, 60 and 100 is labeled Fi2^m ■ ■ ■ Fioo^m- SI is 
the spectral index a determined between 2.7 GHz and 5 GHz. z is the cosmological 
redshift. The distance corrected radio flux flux is given by F^qhI.- The IR data has 
been taken in the 1980s. 
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X-ray source 


Radio source 


Fx 


HR I 


HR II 


F^GHz 


SI 


z 




(2000 coord) 


(1950 coord) 


[cnt / s] 






[Jy] 








IRXS J172320.5+341756 


S4 1721+343 


1.0160 


0.20 


0.07 


0.934 


-0.43 


0.206 


5.2 


IRXS J225358.0+160855 


IJy 2251+158 
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0.96 


0.26 
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0.62 


0.19 
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0.3253 


-0.18 


-0.05 


10.81 


0.54 


0.594 
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Table B.2 

ROSAT sources identified with strong radio blazars. The mean X-ray flux as mea- 
sured by ROSAT in 1990/1991 (6 month of observation time) is given by Fx- The 
hardness ratios HR I and HR II display the ratio between low and high energy X- 
rays: HR = {uhe — ^le) / {nnE + nLs)- For HR I the low energy range is 0.1 — 0.4 
keV and the high energy range is 0.5 — 2.0 keV. For HR II the intervals are 0.5 — 0.9 
keV and 0.9 — 2.0 keV. SI is the spectral index a determined between 2.7 GHz and 
5 GHz. F^GHz stands for the radio flux at 5 GHz. SI is the spectral index a deter- 
mined between 2.7 GHz and 5 GHz with F„ oc The distance corrected radio flux 
flux is given by F^q^^. 
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HEAO-A source 

(1950 coord) 


radio source 

(1950 coord) 


Fx 
[cnt / s] 


F^Ghz 

[Jy] 


SI 


z 


r?z=0.l 

[Jy] 


1H0717+714 


IJy 0716+714 


0.0030 


1.12 


0.22 


0.30 


14.7 


1H1154+294 


S3 1156+29 


0.0029 


0.89 


-0.44 


0.73 


167 


1H1922+746 


IJy 1928+73 


0.0027 


3.34 


-0.01 


0.36 


73.1 



Table B.3 

HEAO-A sources identified with strong radio blazars: The mean flux measured by 
HEAO-A is given by Fx, and F^qhz stands for the radio flux at 5 GHz. AH data 
have been taken within an 6 month interval in 1977-1978, though the observation 
times of individual sources arc shorter. SI is the spectral index a determined between 
2.7 GHz and 5 GHz with ex v°'. The cosmological redshift is given by z and the 
distance corrected radio flux flux is given by F^q^^. 
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Table B.4 

3EG sources identified with strong radio blazars: F^ghz is the radio flux at 5 GHz, 
SI is the spectral index a determined between 2.7 GHz and 5 GHz. The maximum 
and the mean flux as measured by EGRET in 1991-1995 are listed under Fq^y resp. 
FqIy^. z is the cosmological redshift. The EGRET catalog uses the reference system 
of the yea 2000. 
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Table B.5 

Unidentified and weakly identified extragalactic EGRET sources, possibly blazars. 
SI is the spectral index a determined between 2.7 GHz and 5 GHz. F^ghz is the 
radio flux at 5 GHz, SI is the spectral index a determined between 2.7 GHz and 5 
GHz. The maximum and the mean fiux as measured by EGRET in 1991-1995 are 
listed under Fg^ resp. Fg|^°. 

34 



X-ray source 


other name 


FreV 


5 GHz flux 


z 






[crab flux] 


[Jy] 




lES 1101+38.4 


Mrk 421 


0.04-7.4 


0.72 


0.0300 


lES 1652+39.8 


Mrk 501 


0.33-6 


1.38 


0.0336 


H 1426+428 




0.03-0.08 


0.04 


0.1290 


lES 1959+650 




0.05-2.20 


0.25 


0.0470 


lES 2344+514 




0.03 


0.23 


0.0440 



Table B.6 

TeV blazars with 6 > 10°, The TeV fluxes are hsted in units of the flux of the Crab 
nebula, dJ^/dE = (2.79 ± 0.02 ± 0.5) • io-7(^/rey)-2-59±o.03±o.05^ p^^^^ ^^^^^^ 

for the radio flux at 5 GHz, while the cosmological redshift is given by z. 
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Table B.7 

GPS and CSS sources fulfilling our selection rules. In Sample, G means GPS and C 
means CSS. Mag is the optical magnitude. The cosmological redshift is given by z, 
while F^GHz stands for the radio flux at 5 GHz. 9 is the angular size of the source. 

is the turnover frequency and P^ghz stands for log of power at 5 GHz (assuming 
isotropic emission). The linear size of the sources is listed in kpc. The data were 
compiled from previous catalogs in 1998. 

36 



Source 


178MHz 


CT 
Oi 


z 



"D A 

K,A 


z 














z74.L) 


1050. 


-0.76 


12.6667 


12.4667 


0.004 


o4.U 


0I.3 


-0.78 


A'K 1 

41.3167 


3.26667 


c\ r\i T 

0.017 




5d.z 


-0.75 


24.85 


21.35 


0.101 


olU.L) 


55.1 


-0.92 


26.2 


15.0333 


0.054 


660.0 


4d.9 


-1.19 


39.65 


16.4333 


0.029 


4o5.U 


37.8 


-0.75 


26.75 


23.5833 


0.029 


^L- 06.1 


ZD.O 


-0.64 


41.6667 


00000 

6.16666 


0.025 


2d4.U 


ZD.O 


-0.75 


1 n 0000 

19.8833 


11.7 


0.020 


6Lj do.U 


Z4.D 


-0.62 


4:2.7 b 


2.66666 


0.021 




z3.9 


-0.59 


17.15 


18.6 


0.017 


ZYZ.l 


19.4 


-0.60 


13.15 


12.3667 


0.003 


0/^ 000 n 


ion 
18.9 


-0.85 


39.1 


13.6 


0.246 


olo.U 


17.8 


-0.72 


26.3 


15.1833 


r\ 1 no 
0.108 


oU ol.U 


16.8 


-0.57 


00 1000 

32.1333 


1.06667 


0.016 


3C zo.U 


16.3 


-1.06 


. •) / ■ 1000 

26.1333 


n 000000 

0.883333 


0.195 


3C 442. U 


16.1 


-0.96 


10 r 

13.5833 


22.2 


0.026 


3C 305.0 


15.7 


-0.85 


63.4667 


14.8 


0.041 


0/^ 001 r» 

zdl.O 


14.0 


-0.28 


69.9167 


9.85 


C\ C\C\C\ 

0.000 


z9d.O 


13.0 


-0.67 


11.0333 


1 A 0000 

14.2333 


0.023 


z9o.O 


lz.7 


-0.45 


31.6833 


1 0000 

13.8333 


0.045 




IZ.Z 


n 77 


ID. ZOOO 






U.UoZ 


3C 449.0 


11.5 


-0.58 


39.1 


22.4833 


0.017 


3C 346.0 


10.9 


-0.52 


17.35 


16.6833 


0.161 


3C 314.1 


10.6 


-0.95 


70.95 


15.1667 


0.119 



Table B.8 

FR-I radio galaxies: F-nsMHz i« the mean radio flux at 178 MHz, SI is the spectral 
index a determined between 178 MHz and 750 MHz, RA and 6 stand for right 
ascension and declination, respectively in the reference system of 1950. All data 
were collected before 1985. Note, in the original 3CRR catalog, the spectral index 
is defined with an opposite sign. The cosmological redshift is given by z. 
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Table B.9 

FR-II radio galaxies: FnsMHz is the mean radio flux at 178 MHz, SI is the spectral 
index a determined between 178 MHz and 750 MHz RA and 6 stand for right 
ascension and declination, respectvely in the reference system of 1950. All data 
were collected before 1985. Note, in the original 3CRR catalog, the spectral index 
is defined with an opposite sign. The cosmological redshift is given by z. 
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Table B.IO 

Radio- weak quasars from the BQS sorted according to IRAS mean flux at A = 60/xm, 
Feo/im- The cosmological redshift is given by z. The logarithm of the energy flux 
density at 10^^-^ is given by log 2^/14.5. For the BQS, data were taken until 1983, 
IRAS data have been taken in the 1980s. 
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Radio source 


other name 


Samples (Position) 


IJy 0316+41 


NGC 1275 


IR blazars (1), keV blazars (ROSAT (1) and HEAO-A (1)) 


IJy 1652+398 


Mrk 501 


keV blazars (ROSAT (2) and HEAO-A (2)) 


IJy 1222+13 


M 84 


IR blazars (6) 


IJy 0055+30 


NGC 315 


IR blazars (11) 


IJy 1807+698 




keV blazars (HEAO-A (31)) 



Tabic B.ll 

Nearby intrinsically weak sources. Due to their proximity they nevertheless may 
contribute significantly to the neutrino flux. We propose to analyze them as indi- 
vidual sources. 
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